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摘要 :目的 心 提出 一 种 数据 元 余 信 息 引导 的 低 剂 量 心肌 灌注 CT 成 像 方法 。 方 法 考虑 到 心肌 灌注 CT 图 像 帧 内 含有 丰富 的 结 
构 宛 余 性 , 且 帧 间 具 有 高 度 的 相似 性 ,本 文 提出 基于 非 局 部 均值 滤波 (NLM) 和 全 变 分 (TV) 混 合 框 架 的 惩罚 加 权 最 小 二 乘 
(PWLS) 图 像 重建 模型 ,简称 为 PWLS-aviNLM-TV。 该 模型 利用 了 帧 间 结 构 相 似 性 和 帧 内 数据 宛 余 性 ,能 有 效 消除 重建 图 像 中 
的 噪声 和 伪 影 ,提高 灌注 序列 图 像 帧 内 空间 分 辩 率 与 帧 间 时 间 分 辨 率 。 结 果 PWLS-aviNLM-TV 相 比 PWLS-TV 和 PWLS- 
aviNLM 能 更 好 地 去 除 心肌 灌注 图 像 中 的 噪声 和 伪 影 ,同时 较 好 保持 图 像 边缘 和 细节 信息 ,进而 有 效 区 分 缺 血 心肌 与 正常 心 
肌 。 结论 数据 元 余 信 息 引导 的 重建 算法 可 有 效 改善 低 剂 量 心肌 灌注 CT 成 像 质量 , 更 好 地 为 临床 影像 诊断 服务 。 
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Abstract: Objective In the clinic, myocardial perfusion computed tomography (MPCT) imaging is commonly used to detect 
and assess myocardial ischemia quantitatively. However, repeated scanning on the myocardial region in the cine mode will 
increase the radiation dose for patients. With lowering radiation dose, the quality of images are degraded by noise induced 
artifact, which hampers the diagnostic accuracy. Therefore, in this paper, we propose a redundancy information induced 
iterative reconstruction framework for high quality MPCT images at the case of low dose. Methods MPCT images have 
redundant structural information within frames and highly similarity between adjacent frames. Inspired by the two properties, 
in this work we propose a penalized weighted least-squares (PWLS) model incorporating NLM and TV based hybrid 
constraints, which is referred to as PWLS-aviNLM-TV for simplicity. The proposed algorithm can effectively eliminate noise 
and artifacts by taking into account the similarity between adjacent frames and redundancy information within frames, which 
also can improve spatial resolution within frames and maintain temporal resolution. Results The experimental results on the 
4D extended cardiac-torso (XCAT) phantom and preclinical porcine dataset demonstrates that the PWLS-aviNLM-TV 
algorithm obtains better performance in terms of noise reduction and artifacts suppression than the PWLS-TV and PWLS- 
aviNLM algorithm. Moreover, the proposed algorithm can preserve the edges and detail information thereby efficiently 
differentiate ischemia from myocardium. Conclusion The present redundancy information induced reconstruction algorithm 
can reconstruct high-quality images from low-dose MPCT for better clinical imaging diagnosis. 

Keywords: low-dose; myocardial perfusion computed tomography; non-local means; total variation 


验 室 , 广 东 广州 510515 PRR 


近年 来 ,心血 管 疾病 的 发 病 率 呈 持 续 上 升 趋势 ,对 


冠状 动脉 疾病 (CAD) 患 者 进 和 早期 精准 的 诊断 具有 重 
要 的 临床 意义 2。 心肌 灌注 CT 作为 一 种 非 侵 入 式 的 


CAD 功能 成 像 检 查 手 段 , 不 仅 可 以 提供 冠状 动脉 狭 罕 


的 解剖 学 信息 ,还 能 够 获取 冠状 动脉 的 功能 信息 ,在 临 
床上 受到 越 来 越 多 的 关注 5。 根 据 扫 描 协 议 ,心肌 灌注 
CT 对 心 胜 区 域 进行 连续 扫描 和 导 到 时 间 序 列 图 像 , 并 由 
此 计算 出 血 流 动力 学 参数 图 ,为 诊断 心肌 缺 血 区 域 提供 
可 靠 的 影像 诊断 依据 5"。 然 而 ,由 于 心肌 灌注 CT 动态 
扫描 方式 ,导致 病人 接受 的 辐射 章 Ra 而 过 量 
的 XX 出 线 照 射 会 诱发 癌症 等 遗传 性 疾病 。 | 临床 上 党 
A E 
剂量 ,然而 剂量 的 降低 导致 重建 图 像 噪声 和 伪 影 增加 ， 
图 像 质量 退化 ,影响 临床 诊断 的 精确 性 7" “。 因 此 ,如 何 
在 降低 辐射 剂量 的 条 件 下 ,得 到 满足 临床 诊断 需求 的 心 
肌 灌 注 图 像 已 成 为 CT 成 像 领 域 的 热点 研究 课题 之 一 。 


当前 ， 
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近年 来 ,国内 外 大 量 研究 学 者 针对 低 剂量 心肌 灌注 
CT 成 像 进行 了 相关 研究 。 其 中 ,主流 的 研究 方法 大 
致 归 为 两 大 类 :第 1 类 是 图 像 重建 算法 ,该 类 方法 针对 
投影 数据 的 噪声 特性 进行 统计 建 模 ,构建 图 像 迭 代 重建 
模型 实现 图 像 降 品 ,从 而 获得 优质 的 心肌 灌注 CT 图 
像 。Lin 等 ”提出 一 种 无 分 割 的 多 能 MPCT 重 建 算法 ， 
该 算法 结合 了 物质 衰减 系数 ,XX 射线 源 光 谱 和 局 部 组 织 
类 型 的 先 验 信息 ,更 准确 的 重建 出 “无 伪 影 ”的 灌注 网 
像 ;Bian 等 将 心肌 灌注 CT 图 像 空间 以 及 时 间 结 构 的 
稀 玻 特性 引入 到 迭代 重建 过 程 中 ,在 噪声 伪 影 抑制 和 边 
缘 保持 方面 达到 良好 效果 ;Tao 等 “提出 基于 全 变 分 
(TV) 的 统计 迭代 算法 抑制 低 剂 量 心肌 灌注 CT 图 像 的 中 
的 噪声 和 伪 影 ,第 2 类 是 图 像 去 卷 积 算法 ,通过 对 低 剂量 
CT 序列 图 像 进行 去 卷 积 ,从 而 获得 优质 的 血 流 动力 学 参 
数 图 。Zeng 等 人 构造 的 基于 结构 张 量 TV 正则 化 用 于 灌 
注 去 卷 积 算法 ,由 于 TV 算法 关于 图 像 分 段 光 滑 的 假设 ， 
导致 所 重建 的 图 像 中 不 可 避免 的 引入 阶梯 伪 影 *。 

上 述 现 有 算法 均 忽 视 心 肌 灌注 图 像 帧 间 结 构 相 似 
性 和 帧 内 数据 见 余 性 。 鉴 于 此 ,本 文 提 出 数据 元 余 信 息 
引导 的 统计 迭代 算法 ,充分 考虑 灌注 图 像 空间 维度 上 的 
结构 元 余 性 和 时 间 维 度 上 的 信息 相似 性 ,将 心肌 灌注 中 
前 后 相 邻 帧 图 像 的 均值 作为 NLM 算 法 的 参考 图 像 
(aviNLM) ,构造 基于 aviNLM 和 TV 混合 框架 的 惩罚 加 
权 最 小 二 乘 (PWLS)PWLS 模 型 。 采 用 XCAT 体 模仿 真 
数据 和 猪 心肌 灌注 数据 进行 实验 ,与 传统 单 PWLS-TV 
和 PWLS-aviNLM 算 法 相 比 ,PWLS-aviNLM-TV 算 法 
更 能 有 效 地 抑制 低 剂 量 下 的 噪声 和 伪 影 ,同时 较 好 的 保 
持 图 像 结 构 信息 ,进而 有 效 区 分 缺 血 心肌 。 


1 方法 
1.1 惩罚 加 权 最 小 二 乘 CT 重建 模型 

基于 惩罚 加 权 最 小 二 乘 的 低 剂 量 心肌 灌注 CT 统 
计 迭 代 重 建 模型 可 表达 如 下 : 

u*= arg min{(y- Hu) A(y- Hu) +BR(u)} (1) 

其 中 ,为 待 重建 的 心肌 灌注 图 像 ; y 为 经 过 系统 
校准 和 对 数 变换 后 获得 的 投影 数据 ; 五 为 CT 成 像 系统 
矩阵; R(u) 为 正则 化 项 ; 8 为 正则 化 参数 ,其 作用 为 控 
制 R(w) 的 正则 化 力度 。 和 为 权 值 矩 阵 , 其 对 角 元 素 为 
投影 数据 噪声 方差 的 倒数 ,其 余 元 素 为 零 : 

A=diag[ Lt 人 (2) 

o; 

其 中 , c 为 心肌 灌注 投影 数据 中 第 ;个 像素 的 噪声 
方差 ,根据 我 们 之 前 的 研究 ” ,投影 数据 的 噪声 方差 可 
由 如 下 公式 求 得 : 

c? = Looi + TOPP) - 1.25)] (3) 


Hh, 1, Fea ASE KIRE, 元 为 投影 数据 
在 第 i 个 探测 单元 上 的 均值 , of 为 电子 噪声 方差 。 
1.2 双 正 则 化 项 先 验 模型 

为 探索 并 利用 心肌 灌注 图 像 帧 内 宛 余 性 与 帧 间 相 
似 性 特征 ,本 文 将 多 帧 图 像 的 平均 图 像 作 为 先 验 图 像 ， 
并 将 之 纳入 先 验 图 像 引 导 的 非 局 部 均值 正则 化 模型 ,以 
描述 灌注 数据 帧 内 匈 余 性 。 除 此 之 外 ,针对 帧 间 的 高 度 相 
似 性 ,本 文 利 用 全 变 分 正则 化 描述 了 图 像 的 稀 琉 特点 ， 
从 而 构建 “aviNLM-TV” 双 正则 化 模型 ,其 公式 表达 为 : 

R nom- (U) =o plu -aviNLM (u)) 

+(1-B)TVW) 

其 中 , 6,0) WARR, HIERHEEN PTA, ASIC 
以 2 范 数 作为 势 函 数 , 即 P=2 。 pe [0,1] 为 权重 因子 。 

在 心肌 灌注 CT 动态 成 像 中 ,除了 因 造 影 剂 而 增强 
的 血 流动 态 信息 ,序列 网 像 的 前 后 帧 的 结构 具有 高 度 的 
相似 性 ;此 外 ,其 前 后 帧 之 间 的 噪声 是 独立 且 随 机 分 布 
的 。 鉴 于 此 ,本 文 将 当前 图 像 的 前 后 帧 进行 平均 , 取 其 
均值 图 像 作 为 先 验 图 像 并 将 之 纳入 先 验 引 导 的 非 局 部 
均值 正则 化 模型 。 式 (3) 中 , aviNLM (-) 为 均值 图 像 
引导 的 非 局 部 均值 先 验 项 ”2 ,其 表达 式 如 下 : 

aviNLM (u(i)) = ¥o(i,/)u,,(/) (5) 


其 中 , N, 为 搜索 窗 , ulj) 为 均值 参考 图 像 。 
co(i.j) 为 权重 ,其 值 由 低 剂 量 图 像 中 以 像素 i 为 中 心 的 
图 像 块 与 先 验 图 像 中 以 像素 j 为 中 心 的 图 像 块 之 间 的 
相似 性 决定 ,其 计算 公式 为 : 


z ; -C ow; ? 
a(n) -Cu L n 


(4) 


o(i,j) = zy exp 


he 


Z = >》 exp (7) 


JEN; 


P 


$ u(n,)- Cusl), 


Hep n, 和 分别 表 示 为 以 i 和 为 中 心 的 图 像 
块 。|||， 为 两 个 图 像 块 之 间 经 过 高 斯 加 权 的 欧 氏 距 
离 , a 为 高 斯 国 数 的 标准 差 。 式 (6) 中 , 有 h 为 可 调 参 数 ， 
控制 aviNLM 先 验 项 的 平滑 力度 。C 为 局 部 修正 系数 ， 
用 于 修正 待 求 图 像 与 先 验 图 像 之 间 的 局 部 灰 度 水 平 差 
异 , 其 计算 公式 为 : 
C(u(n,),u,,(7,)) = 

E(u,,(n,)) 

1 o th erwise 

其 中 E(-) ARBRE, BUERT y 由 图 像 
块 灰 度 值 的 方差 所 决定 。 

此 外 , TV(-) 表 示 全 变 分 正则 化 先 验 ,其 计算 公式 为 : 


E(u(n) -u,,(n,)) 27 (8) 


H 


TV (u) = > Mu. - ipaa). + (u,, - tics) (9) 
s File 代表 待 重建 心肌 灌注 图 像 & 的 位 置 坐标 。 
1.3 代价 函数 与 优化 求解 
本 文 在 PWLS 框 架 下 引入 了 “aviNLM-TV” 双 正则 
化 先 验 ,代价 函数 具体 表达 如 下 : 


u*=arg min 
wz0 


(y- Hu)'A(y- Hu) +B 
| u- aviNLM (Co + (1 -8)TV (u) 
本 文 利 用 改进 的 梯度 下 降 法 算法 ”对 于 代价 函数 
式 (10) 进 行 迭 代 优 化 求解 。 其 图 像 更 新 公式 为 : 
u” =u" - al (H(X "(Hu -y))) 
-pb Rn rv (u ") 
其 中 a! 表示 梯度 步 长 , a”'= 


(10) 


(11) 
GG 
(HG) (X"(HG)) 
GA 太 (2(B -让 。7z 为 当前 迭代 步 数 , Rw w(x) 
表示 双 正 则 化 项 的 梯度 : 
R om) =p 2, u" - aviNIM (u’)| + 
E pSV (u') 


其 中 如 7F(w) 为 TV 项 的 一 阶 偏 导 , 其 数学 表达 


(12) 


图 1 心肌 灌注 CT 图 像 
Fig.1 Myocardial perfusion CT images. A: XCAT myocardial perfusion image; B: The preclinical porcine 
myocardial perfusion image. 


1.6 动物 数据 

采用 GE 公司 Discovery CT750 64 层 CT 对 猪 进 行 
心肌 灌注 扫描 。 其 扫描 参数 为 : 管 电压 120 kVp, 层 厚 
为 5 mm, 射 源 到 扫描 物体 和 探测 器 的 距离 分 别 为 
538.52 mm 和 946.746 mm, 探 测 器 数目 为 888, 机 架 旋 
转 1 周 为 0.4 s, 旋 转角 度 为 360", 采 样 角度 数 为 984。 同 
上 , 低 剂 量 心肌 灌注 数据 由 高 剂量 数据 仿真 所 得 2”。 


2u,, 一 U,. 1,t 一 U1 


Ou Ie - ial y A i heal rA 
u, 


Usi “Ast 
i 5 5 (13) 


(uu a u,,) 1 Car “Usa. i) +t 


Usti Us 


ION - u) + (ts - tam) +T 

1.4 实验 数据 的 获取 

为 了 验证 所 提 方 法 在 低 剂量 心肌 灌注 CT 成 像 中 
的 有 效 性 和 准确 性 ,本 文 利 用 XCAT 体 模仿 真 数据 二 2 
和 动物 数据 进行 实验 ,并 对 实验 结果 进行 了 定性 分 析 和 
定量 分 析 。 
1.5 XCAT 体 模仿 真 数据 

我 们 首先 根据 血 流动 力学 模型 ** 仿 真 生成 XCAT 
动态 心肌 灌注 序列 图 像 数 据 , 其 中 第 15 帧 心肌 灌注 图 
像 如 图 1A 所 示 。 然 后 通过 正 癌 投影 生成 无 噪声 的 投影 
数据 ,具体 成 像 参数 为 : 层 厚 为 1 mm, 图 像 大 小 为 512x 
512, 射 源 到 探测 需 和 旋转 中 心 的 距离 分 别 为 1040 mm 
和 540 mm, 探 测 器 数目 为 672 ,采样 角度 数 为 1160 ,并 
根据 参考 文献 [27] 中 的 方法 生成 对 应 的 低 剂量 心肌 灌 
注 数据 。 


2 结果 
2.1 XCAT 体 模 数 据 实验 结果 

图 2 给 出 了 不 同方 法 下 第 2 帧 .第 15 帧 和 第 26 帧 
XCAT 心 肌 灌注 图 像 的 重建 结果 。 表 1 列 出 了 图 2 中 不 
同方 法 重建 结果 的 峰值 信 噪 比 (PSNR)。 

图 3 绘 出 了 对 应 于 PWLS-TV,PWLS-aviNLM， 
PWLS-aviNLM-TV 算 法 重建 结果 的 剖面 线 ,剖面 线 的 
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图 2 体 模 真 值 图 像 和 低 剂 量 下 不 同方 法 重建 结果 


PWLS-aviNLM-TV 


Fig.2 True phantom and low-dose images reconstructed by different algorithms. 


表 1 不 同 重建 算法 的 PSNR(dB) 对 比 
Tab.1 PSNR (dB) of different reconstruct algorithms 
Methods FBP PWLS-TV PWLS-aviNLM PWLS-aviNLM-TV 


2 16.65 29.88 26.69 37.44 
15" 16.36 30.43 26.73 37.94 
26" 17.13 30.75 26.84 38.36 


位 置 如 图 3 中 蓝 色 水 平 线 所 示 。 图 4 所 示 为 从 不 同 算法 
重建 的 低 剂 量 序列 图 像 中 计算 得 到 的 心肌 血 流量 
(MBF). 
2.2 动物 数据 实验 结果 

图 5 显示 了 猪 的 第 6 帧 ,第 10 帧 ,第 14 帧 高 剂量 心 
肌 灌注 图 像 和 不 同 算法 的 低 剂量 心肌 灌注 重建 图 像 。 
图 6 所 示 为 猪 的 高 剂量 MBF 参 数 图 以 及 从 不 同 算法 重 
建 的 低 剂量 序列 图 像 计算 得 到 的 MBF 参 数 图 像 。 图 7 
给 出 各 算法 重建 的 低 剂量 图 像 质 量 检测 系数 (UQI)。 


3 讨论 

常规 迭代 重建 算法 的 正则 化 项 设计 往往 只 涉及 图 
像 单一 框架 下 的 先 验 信息 ,例如 , 先 验 图 像 引 导 NLM 滤 
波 的 统计 迭代 重建 算法 ;利用 正常 剂量 图 像 中 的 先 验 信 
息 进行 低 剂 量 重建 ,从 而 提高 重建 图 像 的 质量 ;基于 图 
像 时 间 序 列 先 验 的 正则 化 项 来 重建 心肌 灌注 图 像 3。 
然而 针对 低 剂量 心肌 灌注 CT 成 像 ,本 文 考虑 了 帧 内 结 
构 相 似 性 和 帧 间 相 关 性 ,利用 数据 元 余 信 息 , 提 出 了 一 
种 基于 TV 和 aviNLM 混 合 框架 下 迭代 重建 算法 ,该 算 
法 充分 利用 图 像 自身 分 段 光滑 的 特性 和 心肌 灌注 序列 
图 像 中 不 同 帧 内 结构 相似 性 和 帧 间 相关 性 ,通过 平均 前 
后 帧 的 心肌 灌注 图 像 ,将 帧 间 的 随机 噪声 减少 ,进而 得 
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图 3 不 同 算法 重建 结果 的 剖面 线 比较 


Fig.3 Profile comparison between different algorithms. 


到 较 优 的 参考 图 像 指 导 NLM 去 噪 , 并 将 改进 后 的 
aviNLM 结 合 TV 模型 作为 正则 化 项 融入 PWLS 的 设计 
框架 中 ,构建 出 PWLS-aviNLM-TV 重建 模型 。 对 该 模 
型 的 求解 优化 ,我 们 采用 梯度 下 降 法 算法 来 实现 。 同 
时 ,需要 指出 的 是 ,和 友 代 重 建 过程 中 各 个 参数 的 选择 是 
一 个 公开 的 难 优化 的 问题 ,有 待 今后 更 进一步 的 研究 。 
为 了 评价 PWLS-aviNLM-TV 算 法 的 有 效 性 和 准确 性 ， 
本 文采 用 XCAT 体 模仿 真 数 据 和 猪 心肌 灌注 数据 进行 
了 验证 。 

在 XCAT 体 模仿 真 数据 实验 中 ,图 2 的 FBP 重 建 结 
果 可 以 看 出 ,由 于 剂量 的 降低 ,重建 图 像 存在 大 量 的 伪 
影 和 噪声 ,难以 辨别 出 缺 血 心肌 区 域 ,而 PWLS-TV 和 
PWLS-aviNLM 算 法 重建 结果 中 的 伪 影 和 噪声 大 大 降 
低 , 但 是 仍 存在 阶梯 伪 影 或 者 边缘 模糊 的 现象 。 本 文 提 
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图 4 真 值 和 不 同 算法 重建 低 剂量 图 像 的 MBF 血 流动 力学 参数 图 


Fig.4 MBF calculated from the true phantom and low-dose images reconstructed by different algorithms. 
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图 5 高 剂量 心肌 灌注 图 像 和 不 同 算法 重建 的 低 剂 量 心肌 灌注 图 像 


Fig.5 High-dose and low-dose myocardial perfusion images reconstructed by different algorithms. 
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图 6 高 剂量 MBF 参数 图 和 不 同 算法 重建 下 低 剂 量 图 像 的 MBF 参 数 图 


Fig.6 MBF maps calculated from a high-dose image and simulated low-dose images reconstructed by different algorithms. 


出 的 双 正 则 化 方法 的 重建 结果 能 有 效 保持 边缘 结构 信 ”和 较 单一 正则 化 TV 算法 ,单一 aviNLM 算 法 和 双 正 则 化 
息 , 并 抑制 伪 影 和 噪声 。 从 表 1 定 量 分 析 上 可 以 看 出 ， ”算法 对 边缘 结构 的 保持 效果 ,图 3 给 出 重建 结果 剖面 线 
本 文 提 出 的 算法 重建 图 像 的 PSNR 高 于 其 它 三 类 算法 ”的 比较 ,其 结果 显示 , 相 比 TV 和 aviNLM 单 一 正则 化 算 
重建 图 像 ,能 够 有 效 的 抑制 了 图 像 噪声 。 为 了 进一步 比 。 法 ,PWLS-aviNLM-TV 算 法 与 真 值 的 剖面 线 最 为 接近 ， 
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图 7 高 剂量 心肌 灌注 图 像 和 不 同 算法 重建 的 低 剂量 心肌 


灌注 图 像 的 UQI 

Fig.7 UQI measurements on a high-dose myocardial 
perfusion image and simulated low-dose myocardial 
perfusion images reconstructed by different algorithms. 


即 PWLS-aviNLM-TV 算 法 能 有 较 好 的 重建 出 边缘 结 
构 。 为 了 评估 重建 算法 的 准确 性 ,我 们 采用 经 典 的 奇异 
值 分 解 方法 汪 计 算 了 心肌 血 流动 力学 参数 图 ,如 图 4 
所 示 , 可 以 看 出 FBP 算 法 重建 的 图 像 存在 大 量 的 噪声 和 
伪 影 ,以 至 于 无 法 辨别 缺 血 心肌 ,而 PWLS-TV 和 
PWLS-aviNLM 算 法 的 结果 稍 有 改善 ,但 与 真 值 图 像 还 
有 一 定 的 差距 , 相 比 之 下 ,本 文 PWLS-aviNLM-TV 的 
方法 在 抑制 噪声 的 同时 能 保证 较为 精确 的 MBF 参 数 图 
估计 。 

在 猪 心肌 灌注 数据 实验 中 ,图 5 的 结果 可 观察 到 ， 
本 文 提出 的 PWLS-aviNLM-TV 算 法 不 仅 可 以 有 效 抑 
制图 像 的 噪声 和 伪 影 ,同时 较 好 保持 了 具有 诊断 意义 的 
结构 信息 。 从 MBF 参 数 图 定性 分 析 结 果 中 可 以 看 出 ， 
FBP 算 法 重建 的 图 像 中 存在 大 量 噪声 和 伪 影 ,导致 其 
MBF 参数 图 无 法 辨别 心肌 缺 血 区 域 ,另外 ,PWLS-TV 
和 PWLS-aviNLM 算 法 重建 的 图 像 一 定 程度 上 可 辨别 
心肌 缺 血 区 域 ,但 是 在 图 6 白色 虚线 框 的 缺 血 心肌 区 域 
仍 存在 部 分 噪声 ,上 且 边 缘 细节 有 所 退化 , 相 比 之 下 ， 
PWLS-aviNLM-TV 算 法 重建 的 图 像 不 仅 抑制 了 噪声 和 
伪 影 ,而 且 边 缘 细节 与 高 剂量 的 图 像 具 有 较 高 的 一 臻 
性 。 通 过 图 7 定量 结果 显示 ,不同 时间 帧 下 ,UQI2" 得 分 
从 低 到 高 依次 为 FBP,PWLS-TV,PWLS-aviNLM 和 
PWLS-aviNLM-TV 算 法 ,这 进一步 表明 ,本 文 提出 的 
PWLS-aviNLM-TV 算法 重建 的 低 剂量 图 像 与 高 剂量 
图 像 在 结构 和 细节 上 更 为 接近 。 

综 上 所 述 , PWLS-aviNLM-TV 算 法 结合 了 NLM 
和 TV 模型 特性 ,能 更 有 效 地 抑制 低 剂量 心肌 灌注 图 像 
中 的 伪 影 和 噪声 ,提高 灌注 序列 网 像 帧 内 空间 分 辩 率 与 
帧 间 时 间 分 辩 率 ,上 且 从 其 MBF 人 参数 图 中 可 更 清晰 判断 


出 缺 血 区 域 的 位 置 和 大 小 。 
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